An entanglement-area law for general bosonic harmonic lattice systems 
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We demonstrate that the entropy of entanglement and the distillable entanglement of regions with respect 
to the rest of a general harmonic lattice system in the ground or a thermal state scale at most as the boundary 
area of the region. This area law is rigorously proven to hold true in non-critical harmonic lattice system of 
arbitrary spatial dimension, for general finite-ranged harmonic interactions, regions of arbitrary shape and states 
of nonzero temperature. For nearest-neighbor interactions - corresponding to the Klein-Gordon case - upper and 
lower bounds to the degree of entanglement can be stated explicitly for arbitrarily shaped regions, generalizing 
the findings of [Phys. Rev. Lett. 94, 060503 (2005)]. These higher dimensional analogues of the analysis of 
block entropies in the one-dimensional case show that under general conditions, one can expect an area law 
for the entanglement in non-critical harmonic many-body systems. The proofs make use of methods from 
entanglement theory, as well as of results on matrix functions of block banded matrices. Disordered systems are 
also considered. We moreover construct a class of examples for which the two-point correlation length diverges, 
yet still an area law can be proven to hold. We finally consider the scaling of classical correlations in a classical 
harmonic system and relate it to a quantum lattice system with a modified interaction. We briefly comment on 
a general relationship between criticality and area laws for the entropy of entanglement. 

PACS numbers: 03.67.Mn, 05.50. +q, 05.70.-a 
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I. INTRODUCTION 

Ground states of quantum systems with many constituents 
are typically entangled. In a similar manner as one can iden- 
tify characteristic length scales of correlation functions, quan- 
tum correlations are expected to exhibit some general scal- 
ing behavior, beyond the details of a fine-grained descrip- 
tion. Such characteristic features provide a physical picture 
that goes beyond the specifics of the underlying microscopic 
model. A central question of this type is the following: If one 
distinguishes a certain collection of subsystems, representing 
some spatial region, of a quantum many -body system in a pure 
ground state, the state of this part will typically have a posi- 
tive entropy, reflecting the entanglement between this region 
and the rest of the system fl |3H @, H H 0] . This de- 
gree of entanglement is certainly expected to depend on the 
size and also on the shape of the region. Yet, how does the 
degree of entanglement specifically depend on the size of the 
distinguished region? In particular, does it scale as the volume 
of the interior - which is meant to be the number of degrees 
of freedom of the interior? Or, potentially as the area of the 
boundary, i.e., the number of contact points between the inte- 
rior and the exterior? 

This paper provides a detailed answer to the scaling be- 
havior of the entanglement of regions with their exterior in 
a general setting of harmonic bosonic lattice systems and pro- 
vides a comprehensive treatment of upper bounds on these 
quantities. We find that in arbitrary spatial dimensions the 
degree of entanglement in terms of the von Neumann entropy 
scales asymptotically as the area of the boundary of the distin- 
guished region. This paper significantly extends the findings 
of Ref. |sj on harmonic bosonic lattice systems: There, the 
area-dependence of the geometric entropy has been proven 
for cubic regions in non-critical harmonic lattice systems of 
arbitrary dimension with nearest-neighbor interactions, cor- 



responding to discrete versions of Klein-Gordon fields. In 
this paper we extend our analysis to a general class of finite- 
ranged harmonic interactions and also take regions of arbitrary 
shape into account. For thermal Gibbs states, the entropy of 
a reduction is longer a meaningful measure of entanglement. 
Instead, an area-dependence for an appropriate mixed-state 
entanglement measure, the distillable entanglement, is estab- 
lished. Also, an analogous statement holds for classical cor- 
relations in classical systems. The area-dependence is even 
found in certain cases where one can prove the divergence of 
the two-point correlation length. This demonstrates that this 
previously conjectured dependence between area and entan- 
glement is valid under surprisingly general conditions. 

The presented analysis will make use of methods from the 
quantitative theory of ent ang lement in the context of quantum 
information science 

(minim. 

It has been become clear 
recently that on questions about scaling of entropies and de- 
grees of entanglement - albeit often posed some time ago - 
new light can be shed with such methods fl B |3 l B, E |l a Bl p 

u m Hi m in m m m m eh es ihh ninii. 

In this language, quantum correlations are sharply grasped in 
terms of rates that can be achieved in local physical transfor- 
mations. To assess quantum correlations using novel pow- 
erful tools from quantum information and to relate them to 
information-theoretical quantities constitutes an exciting per- 
spective. 

In the context of quantum field theory, such questions of 
scaling of entropies and entanglement have a long tradition 
under the keyword of geometric entropy. In particular, work 
on the geometric entropy of free Klein-Gordon fields was 
driven in part by the intriguing suggested connectio n 12811 to 
the Bekenstein-Hawking black hole entropy EH [fir In 
seminal works by Bombelli et al. l32ll and Srednicki |33] the 
relation between the entropy and the boundary area of the re- 
gion has been suggested and supplemented with numerical ar- 



guments. This connection has been made more specific using 
a number of different methods. In particular, for half spaces 
in general and intervals in the one-dimensional case, the prob- 
lem has been assessed employing methods from conformal 
field theory, notably 1341 13511 . based on earlier work by Cardy 
and Peschel l36ll . and by Cardy and Calabrese ffl . 

In one-dimensional non-critical chains, one observes a sat- 
uration of the entanglement of a distinguished block, as was 
proven analytically for harmonic chains 0] and was later 
observed for non-critical spin chains numerically 01 and an- 
alytically 171, Il3l. FA . In turn, in critical systems systems, one 
often - but not always - finds a logarithmically diverging en- 
tropy 00 13 13 In case of a model the continuum limit of 
which leads to a conformal field theory, the factor of the log- 
arithmically diverging term is related to the central charge of 
the conformal field theory. The findings of the present paper 
and of the above-mentioned results motivate further questions 
concerning the general area-dependence of the degree of en- 
tanglement, for example in fermionic systems 1241 12511 . In par- 
ticular, the connection between the geometric entropy fulfill- 
ing no area law and the correlated quantum many -body system 
being critical is not fully understood yet. This is particularly 
true for the interesting case of more than one-dimensional 
quantum systems. 

This paper is structured as follows. We start, in Section ITTI 
with a presentation of the major results of the paper. In Sec- 
tion|^we define our notation and recall some basics on har- 
monic lattice systems. Section ffvl provides a general frame- 
work of upper and lower bounds for entanglement measures, 
expressed in terms of the spectrum of the Hamiltonian and 
two-point correlation functions. This analysis is performed 
both for the case of the ground state as well as for mixed 
Gibbs states. Of particular interest are Hamiltonians with 
finite-ranged interactions in Section |V] For such Hamilto- 
nians, we first study the behavior of the two-point correla- 
tion functions, then the entanglement bounds. In Section IVT1 
discussing the Gibbs state case, we determine temperatures 
above which there is no entanglement left. A class of exam- 
ples of Hamiltonians which exhibit a divergent two-point cor- 
relation length in their ground state, but an area-dependence 
of the entanglement is presented in Section Hxl and expressed 
in analytical terms. The specific case of Hamiltonians the 
interaction part of which can be expressed as a square of a 
banded matrix is discussed in Section IVIIII In this case, very 
explicit expressions for entanglement measures can be found. 
We then consider, in Section^ the case of classical correla- 
tions in classical harmonic systems with arbitrary interaction 
structure. Interestingly, this case is related to the quantum case 
for squared interactions in the sense of Section lVIIII Finally, 
we summarize what has been achieved in the present paper, 
and present a number of open questions in this context. 



II. MAIN RESULTS 

Throughout the whole paper we will consider har- 
monic subsystems on a D-dimensional cubic lattice C = 
[1, . . . , n] xD . The system thus has n D canonical degrees of 
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FIG. 1: A two-dimensional lattice system C = [1, n] x with a 
distinguished (shaded) region /, consisting of v(I) degrees of free- 
dom •. Oscillators representing the exterior O — C \ I are shown 
as O, whereas pairs belonging to the surface s(J) of the region are 
marked by lines (C*^). 



freedom. The central question of this work will be: how does 
the degree of entanglement of a distinguished region I C C 
with the rest O = C\ I scale with the size and shape of II 

We define the volume v(I) and surface s(I) of the distin- 
guished region I as 

«(/)= J> «(!)=£ E L 



iei 



where 



D 



ieo jei 
d(i,j)=i 



*<5 - 3s\ 



5=1 



defines the 1-norm distance for vectors i = (ix, . . . , in) € C 
that specify the position of oscillators on the D-dimensional 
lattice. More specifically, s(7) is the number of contact points, 
that is, the number of pairs of sites of I and O that are immedi- 
ately adjacent. Note that the distinguished region I may have 
arbitrary shape and does not have to be contiguous. 

For the pure ground state g, we will study the entropy of 
entanglement of I C C given by 

where S = — tr[p\ogp] is the von Neumann entropy of a 
state p and gj = tro[o\ denotes the reduced state associ- 
ated with the degrees of freedom of the interior I. For the 
pure ground state, this entropy of entanglement is identical to 
both the distillable entanglement and the entanglement cost. 
For pure states, it is indeed the unique asymptotic measure of 
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entanglement. For Gibbs states, in turn, we have to study a 
mixed-state entanglement measure, as the entropy of a sub- 
system no longer meaningfully quantifies the degree of entan- 
glement. We will bound the rate at which maximally entan- 
gled pairs can be distilled, i.e., the distillable entanglement 
E%j(T). Clearly, Egj(0) = E^j for zero temperature. 

We will subsequently suppress the index C for notational 
clarity. We derive the following properties of Ef and Ef (T): 

(I) For D- dimensional harmonic lattice systems, we derive 
general upper and lower bounds to Ef for pure ground 
states and to Ef (T) for Gibbs states with respect to 
a temperature T > 0. These bounds are expressed 
entirely in terms of the potential matrix and stated in 
Eqs. (0, and r TTOt . 

A necessary condition for the following results to hold is that 
the spectral condition number k = A max (l/)/A m i n (F) of the 
coupling matrix V satisfies k < c < oo for some c > 
independent of I and O. 

(II) For nearest-neighbor interactions and the ground state, 
the entropy of entanglement Ef scales as the surface 
area s(I) of I. More specifically, there exist numbers 
C\ , C2 > independent of O and I such that 

ci • s(7) < Ef < c 2 • s(7). 

This is stated in Eqs. M2\ and M5\ . Note that the spe- 
cific case of cubic regions I = [1, . . . , m] xD has al- 
ready been proven in the shorter paper Ref. J^/. 

(III) For general finite-ranged harmonic interactions - 
meaning arbitrary interactions which are strictly zero 
after a finite distance - the entropy of entanglement Ef 
of the ground state scales at most linearly with the sur- 
face area of I: There exists a c > independent from 
O and I such that 

Ef <c-s(I). 

This is expressed in Eq. M2\ . 

(IV) For Gibbs states with respect to a temperature T > 
and general finite-ranged interactions the distillable en- 
tanglement Ef (T) scales at most linearly with the sur- 
face area of I, so 

Ef(T)<c(T)- S (I), 

with c(T) > being independent of O and I, see again 
Eq. M2\ . which applies for any temperature. We also 
determine temperatures above which the entanglement 
is strictly zero, see Eq. (lU. 

(V) We construct a class of Hamiltonians the ground state 
of which exhibit an infinite two-point correlation length, 
for which the entropy of entanglement Ef scales prov- 
ably at most linearly in the boundary area. This is 
stated in Eq. M 91 . 



(VI) For interactions that are specified by potential matrices 
V (see Eq. (Qj that can be written as V = W 2 with W 
corresponding to a finite-ranged interaction, we deter- 
mine the entropy of entanglement. This is made specific 
in particular for D = 1, see Eq. M81 . 

(VII) For classical systems of harmonic oscillators with 
nearest-neighbor interactions prepared in a thermal 
state, the mutual information 

I I (T) = S I (T) + So(T)-S c {T) 

measuring classical correlations scales linearly with 
the surface area of the region. Here, Si, So, and Sc 
are the discrete classical entropies in phase space of 
the interior, the exterior, and the entire system with 
respect to an arbitrary coarse graining as defined in 
Eq. \22\ . for a classical Gibbs state at temperature 
T > 0. The mathematical formulation of the problem 
is intimately related to case (VI). 

The remainder of the paper is now concerned with the de- 
tailed discussion and rigorous proofs of these statements. 



III. HARMONIC LATTICE SYSTEMS - PRELIMINARIES 

We consider quantum systems on Z?-dimensional lattices, 
where each site is associated with a physical system. The 
starting point is the following Hamiltonian which is quadratic 
in the canonical coordinates, 

i i,j 

The matrix V is the potential matrix. Vi t i denotes the potential 
energy contained in the degree of freedom labeled with i. For 
i j the element Vi.j describes the coupling between the two 
oscillators at i and j. We assume that V is a real symmetric 
positive matrix. In this paper, we will consider finite-ranged 
and nearest-neighbor interactions. Furthermore, interactions 
for which the correlation length diverges are considered. 

We will subsequently study properties of both the ground 
state q(0) as well as of Gibbs states 

e -H/T 
Q{T) = tx[e- H / T ] 

corresponding to some nonzero temperature T. Note that we 
have set the Boltzmann constant k = 1. The questions we ask 
are essentially those of the geometric entropy (T = 0), and 
those on the distillable entanglement (T > 0) with respect to 
a distinguished region I of the lattice C. The specific case 
of nearest-neighbor interactions, considered in Refs. Isll^Eol 
l33ll corresponds to the one of a discrete version of free Klein- 
Gordon fields in flat space-time. Here, we consider regions of 
arbitrary shape and allow for more general interactions. 
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A. Phase space, covariance matrices, and ground states 

The system we consider embodies n D canonical degrees of 
freedom, associated with a phase space (R x2 ™ , a), where 
the (2n D ) x (2 n B )-matrix a, 



a :- 



n D l n D 
— 1,,-D 0„D 



specifies the symplectic scalar product, reflecting the canon- 
ical commutation relations between the 2 n D canonical coor- 
dinates o = (xi , . . . , x n o ,pi,..., p n D ) of position and mo- 
mentum. Instead of considering states we will refer to their 
moments, considering that both the ground state as well as the 
Gibbs states are quasi-free (Gaussian) states. The second mo- 
ments can be collected in the (2 n D ) x (2 n D ) real symmetric 
covariance matrix 7 (see, e.g., Ref. fioll ). defined as 

7j -, fe := 2^(tr[ Q o j0k \) = 2^( 0j o k ) e . 

Here, it is assumed that the first moments vanish, which does 
not restrict generality as they can be made to vanish by lo- 
cal unitary displacements that do not affect the entanglement 
properties. 

Following Ref. Hj, we find from symplectic diagonaliza- 
tion that for the ground state at zero temperature 



To 



y-l/2 e yl/2 



i.e., there is no mutual correlation between position and mo- 
mentum, and the two-point vacuum correlation functions are 
given by the entries of V^ 1 ! 2 respectively V -1 / 2 , 



G itJ := (gsl^-lgs) = [V 1/2 }ij, 



Hi 



(gs\piPj\gs) = [V 



1/21. . 



(2) 

(3) 



Here, |gs) denotes the state vector of the ground state. For the 
thermal Gibbs state at finite temperature the covariance matrix 
takes the form 

7 (r) = (y- 1/2 w{T)} ® (v 1/2 w(T)j , 



where we define 



W{T) := 1 + 2(cxp(y 1/2 /T) - 1) 



(4) 



Note that the additional term W(T) is the same for the posi- 
tion as well as for the momentum canonical coordinates. As in 
the zero temperature case position and momentum are not mu- 
tually correlated and V _1 / 2 W(T) {V X / 2 W{T)) are the two- 
point correlation functions of position (momentum), now with 
respect to the Gibbs state. 



B. Measures of entanglement 

The entropy of entanglement can be expressed in terms of 
the symplectic eigenvalues of the covariance matrix corre- 
sponding to a reduction. Let the (2v(I)) x (2u(/))-matrix 



70 1 j denote the covariance matrix associated with the interior 
/; this is the principle submatrix of 70 associated with the de- 
grees of freedom of the interior. The symplectic spectrum of 
7o| 7 is then defined as the spectrum n(Bi) of the matrix 



Br 



To 



1/2 



( icr l/) (To 



1/2 



For the situation at hand we find that 



To 



where we denote by V ±x l 2 \i the principle submatrix of 
V ±x l 2 associated with the interior and we index the entries of 
the submatrix as \V~ 1 ^ 2 ]i,j with vectors i,j g /. Given the 
direct-sum structure of 70 |j and counting doubly degenerate 
eigenvalues only once, we find that the entropy of entangle- 
ment can be evaluated as 



v(I) 



E) 



£[/( 



-/ 



Mi 



(5) 



where f(x) — — x\ogx, and the \ii — ^i(Aj) are the square 
roots of the regular eigenvalues Xi = Xi(Aj) of the v(I) x 
v(I) matrix Aj, 



IM = VMA7), Aj = V- 1 ' 2 



yV2 



This reflects the fact that the entropy of a harmonic system 
is the sum of the entropies of uncoupled degrees of freedom 
after symplectic diagonalization. 

Using vectors i,jEl to label the entries of Aj, we find due 

to V- X I 2 V Y I 2 = 1, 



v- 1/2 \ V 1 ' 2 
\l 



I 

- *j - E [y- 1/2 

keo 



1,3 



yV2 



(6) 



k.j 



This form of Aj = 1 — R hints at an area theorem in the 
following way (see Fig. [2] where we depict the entries of i?): 
if the entries of V ±x l 2 decay fast enough away from the main 
diagonal, i.e., if the correlation functions decay sufficiently 
fast, the main contribution to the entropy comes from oscilla- 
tors inside a layer around the surface of / as for all the others 
the product [V^^Ji^lV 1 / 2 ]*;^ will be very small (here vec- 
tors i, j are inside /, k outside the region /). Thus, the matrix 
R has an effective rank proportional to s(J), and the number 
of symplectic eigenvalues contributing to the sum in Eq. (f5jl 
is approximately proportional to s(I). Much of the remainder 
of the paper aims at putting this intuition on rigorous grounds. 

For mixed states, such as thermal Gibbs states, the von Neu- 
mann entropy no longer represents a meaningful measure of 
the present quantum correlations, and has to be replaced by 
concepts such as the distillable entanglement. The distillable 
entanglement is the rate at which one can asymptotically dis- 
till maximally entangled pairs, using only local quantum oper- 
ations assisted with classical communication. For pure states, 
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Eq. (|4}. The diagonal matrix P, defined as 




FIG. 2: Entries of R in one dimension, D = 1, for a finite- 
ranged coupling matrix V, C = [1, 100] and 7 = [30, 35] U 
[65, 70], yielding s(I) = 4. Bars show Rij, color-encoded sur- 
face depicts log(Ri,j). All units are arbitrary. Note that the entries 
decay exponentially away from the boundary of I. 



it coincides with the entropy of entanglement - the entropy of 
a reduction - giving an operational interpretation to this quan- 
tity noon!. 

From entanglement theory we know that an upper bound 
for the distillable entanglement is provided by the logarithmic 
negativity 13711 . Note that this is the case even in this infinite- 
dimensional context for Gaussian states can be immediately 
verified on the level of second moments and a single-mode 
description. Moreover, the entropy of entanglement indeed 
still has the interpretation of a distillable entanglement, albeit 
the fact that distillation protocols leave the Gaussian setting. 
This is true as long as one includes an appropriate constraint 
to the mean energy in the distillation protocol ll38ll . 

The logarithmic negativity 13711 is defined as 

E?(T) = \\ e (T) r \\ lt 

where ||.||i denotes the trace norm, and g(T) r is the partial 
transpose of g(T) with respect to the split I and O = C \ I. 
Again following Ref. |5], we find after a number of steps 



Ef{T) = log 2 [max {l, A* (Q) } 



(7) 



where Ai(Q) labels the n D eigenvalues of Q, 

Q := Puj-Plu + , 
and matrices ui^ (T) are defined as 

u ± {T) := VF(T)~V ±1/2 , 
which become = V ±x / 2 for zero temperature following 



P. 



Mes = 



$ij died 



1 ifieS 
otherwise, 



is the matrix that implements time reversal in the subsystem 
corresponding to the inner part I, reflecting partial transposi- 
tion g r on the level of states. 



IV. UPPER AND LOWER BOUNDS 

In this section we will derive upper and lower bounds for 
the entropy of entanglement and the distillable entanglement 
of the distinguished region I with respect to the rest of the 
lattice. These bounds only depend on the geometry of the 
problem, i.e., the region I and on properties of V, namely its 
minimal and maximal eigenvalue, which we define as a : = 
Amin(^) and b := A max (V^), respectively, its condition number 



AmaxOO 



and the entries of V^ 1 / 2 respectively the two-point correlation 
functions G and H as in Eqs. (I2I3I . In later sections these 
bounds will be made specific for a wide range of interaction 
matrices V . 



A. Upper bound 

An upper bound for the entropy of entanglement and the 
distillable entanglement is provided by the logarithmic nega- 
tivity as in Eq. (2J. Utilizing this fact, we derive upper bounds 
using Zi-norms 13911 in this section. For brevity and clarity, 
we will present the case for T = and finite temperature in a 
single argument. 

A direct calculation shows that the matrix Q is given by 

Q = - 2Xuj + 

= W{T)- 2 - 21w+, 

introducing the matrix X with entries 



Xi 



~ x.j \ 1 J ieo°jei ; 

Therefore, we can bound the eigenvalues of Q according to 

HQ) < Kin(W(T))- 2 + X i (-2Xuj+) 
< \ m UW(T))- 2 + 2\\i(Xuj + ) 

where we denote by A rn i n (M / (T)) the smallest eigenvalue of 
W(T) which is given by 



A min (^(T))- 1 = A n 



^ 1/2 / T + 1 



; Vb/T ■ I ■ 
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Hence, we can write 



glement as follows, 



<(T)<^>g 2 



iec 



max l,A min (W(T))- 2 + 2 X i (Xtu+) 



< 



ln^) E max { ' A mi „(iy(T))- 2 - 1 + 2 |Ai(Xw+)|j 



i.e., for A min (iy(T))- 2 + 2ma Xi \Xi{Xui+)\ < 1, there is 
no longer any bi-partite entanglement in the system. We will 
later see that there is a temperature T c above which this hap- 
pens. But for now, we use the fact that \ m m(W(T))~ 2 < 1 to 
bound the logarithmic negativity, and relate it to the l\ norm 
\\.\\ h |H of X. We have A max (w+) = Vb/ X min (W(T)) and 
therefore 

&W<^E*Mx» + )\ = ^\\x<* + \U 



iec 



< 



2y/b 



A min (W(T))ln(2) 



\X\\i< 



ln(2) 1 
2Vb 



A min (W(T))ln(2) 



\X\ 



The li norm is defined as the sum of the absolute values of all 
matrix entries. Inserting the definition of X, we find 



A min (^(r))ln(2) 
2Vb 



E?<\\X\\ h = J2 K 



ijec 



2 E| 

ieo 



(8) 



i,3 



for finite temperature, and 



E 



(9) 



for zero temperature. These constitute upper bounds on the 
entropy of entanglement and the distillable entanglement. 
Both only depend in the distinguished region /, the maximum 
eigenvalue of V, and the entries of U -1 / 2 , i.e., the two-point 
correlation function with entries Gij. This will be the start- 
ing point to derive explicit upper bounds for special types of 
interaction matrices V. 



B. Lower bound 

To achieve lower bounds for the zero temperature case, 
we consider the entropy of entanglement directly. Starting 
from the general expression Eq. (|5}, we use the fact that the 
symplectic eigenvalues are never smaller than 1, i.e., that the 
eigenvalues of Ai are contained in the interval [l,a/], with 
a j being the maximal eigenvalue of Aj. Using the pinch- 
ing inequality JH, we find a/ < (A max (U)/A min (U)) 1/2 = 
(6/a) 1 / 2 = k 1 / 2 . Thus we can bound the entropy of entan- 



«C0 «C0 , rx i 



log[aj] Vf v 



Using vectors to label the entries of the v(I) x v(I) 

matrix Aj, we finally arrive at the lower bound 



Ef > - 



2(V«-1) 



yV2 



(10) 



jeo 



This lower bound depends only of the geometry of /, the spec- 
tral condition number n of V and the entries of V ±x l 2 , i.e., the 
two-point correlation functions. 

Eq. l llOt is difficult to evaluate in general but for special 
cases of the interaction matrix V it can nevertheless be made 
specific. In Section fvTTI for example, we give an explicit ex- 
pression for the important case of nearest-neighbor interac- 
tions. We expect Eq. dlOt to be a convenient starting point to 
derive such lower bounds also for more general cases of V. 

A lower bound for the finite temperature case is difficult to 
obtain. Generally speaking, a lower bound to the distillable 
entanglement (two-way distillable entanglement in our case) 
is given by the hashing inequality l4(ill . 



Ef(T) > max{S , (( ?/ ) 



S(q),S(q o )-S(q),0}, 



where gi = tro [g] and go = tr / [p] . Yet, naively applied, this 
inequality will vanish, as generally S(gj) — S(g) < and 
S(go) — S(g) < 0: this can be made intuitively clear from 
the following argument. The above analysis demonstrates that 
both the interior / and the exterior O can be approximately 
disentangled with local unitaries up to a layer of the thick- 
ness of the two-point correlation length. That is, degrees of 
freedom associated with i 6 O for which d(i,j) is suffi- 
ciently large for every j € s(I), can to a very good approxi- 
mation be decoupled and unitarily transformed into a thermal 
Gibbs state. Each such degree of freedom will therefore con- 
tribute a constant number to S(g), such that S(qi) — S(q) < 
for sufficiently large n. Similarly, one can argue to arrive at 
S{qo) ~ S(g) < 0. In order to establish a lower bound to 
the distillable entanglement, however, we may start with any 
protocol involving only local quantum operations and classi- 
cal communication, and apply the hashing inequality on the 
resulting quantum states. This first step can include in partic- 
ular local filterings. A bound linear in the boundary area is 
expected to become feasible if one first applies an appropriate 
unitary both in / and O, and then performs a local filtering in- 
volving degrees of freedom associated with i E I and j € O 
for which d(i, j) = 1. This option will be explored elsewhere. 



V. FINITE-RANGED INTERACTIONS 

In this section we will make the upper bound on the entropy 
of entanglement and the distillable entanglement explicit for 
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symmetric finite-ranged interaction matrices V, i.e., matrices 
for which 



vs. 



0, for d{i,j) > k/2, 



where d(i,j) denotes again the 1-norm distance. Denoting 
as before the maximum and minimum eigenvalue of V as 
a = Amin(V) and b = X max (V) respectively, we require that 
the spectral condition number k = b/abe strictly less than 
infinity independent of C, i.e., independent of n. Note that we 
do not require any further assumptions on the matrix V. 



A. General upper bounds for finite-ranged interactions 




We will make use of a result of Ref. ITHl concerning the 
exponential decay of entries of matrix functions. After gener- 
alizing this result to matrices V with the properties specified 
above it enables us to bound the entries of uj~ as follows (see 
AppendixlAl 



and 



where 



<K a>b q% 



K a ,b = 



n 



d(i,j) 



K-l 
K + l 



2/k 



(ID 



e v/T + 1 „ ; 



K + l 

For zero temperature we then find 



1/2 



K a , b = 



K+l f K+l 



1/2 



This shows that off-diagonal terms of V^ 1 ^ 2 decay exponen- 
tially, cf. Fig. [3] Substituting Eqs. m t° the general result 
in JHJ, we find 



FIG. 3: The entries [V -1 ' 2 ]<,j for a finite-ranged coupling matrix V 
in two dimensions, D = 2, j = i + (xi, X2). Bars show [V~ 1,l2 ]ij 
and the color-encoded surface shows log([V~ 1/ ^ 2 ]ij). The inset de- 
picts the same for xi = X2- All units are arbitrary. Note the expo- 
nential decay away from the main diagonal i — j. 



Thus we finally arrive at the desired upper bound linear in the 
surface area of / for both the entropy of entanglement and the 
distillable entanglement, 



Ef (T) < 



4VbK„ b K 



A min (W(T))ln(2) 



(12) 



which becomes 



< 



4(k+1)2K d , 
ln(2) 



■s(I) 



for zero temperature, i.e., these upper bounds depend solely 
on the maximal and minimal eigenvalue of the interaction ma- 
trix V and the surface area of the region I. This result demon- 
strates that indeed, an area-bound of the degree of entangle- 
ment holds in generality forbosonic harmonic lattice systems. 
This shows that the previously expressed intuition can indeed 
be made rigorous in form of an analytical argument. 



4\/MiT„,, 



iei 
jeo 



where 



jeo iei 
d(i,j)=i 



1=1 



Note that Ni = s(I) coincides with the definition of the 
surface area of I. We find Ni < 2(2l) 2D - 1 s(I) (see Ap- 
pendixO. i.e., we have 

^4^)<2™s(I)jr q lj™-i=:K D , K s(I). 



iei 
jeo 



1=1 



B. Disordered systems 

Notably, the derived results hold also for systems in which 
the coupling coefficients are not identical, but independent re- 
alizations of random variables. If the coefficients VS ,j of the 
real symmetric matrix V are taken from a distribution with a 
carrier [so,xi], such that the interaction is (i) finite-ranged, 
and (ii) the carrier is chosen such that 



A min (V) 



< 00, 



then the same result holds true. This follows immediately 
from the considerations in Appendix [A] where it is not as- 
sumed that the Hamiltonian exhibits translational symmetry. 
Eq. dl 2i is thus valid also for disordered harmonic lattice sys- 
tems. 
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VI. TEMPERATURES ABOVE WHICH THERE IS NO 
MORE ENTANGLEMENT LEFT 



In Section HVAl we found that for X min {W (T))- 2 + 
2maxj |Ai(Xoi + )| < 1 there is no entanglement between the 
regions / and O. For finite-range interactions V we found in 
Section[V] 

max|A i (Xw+)| < XMWiT^Vbn, 

i 

< 2\ min (W(T))- 1 VbK a , b K D:K s(I), 

i.e., we have 

A m i„(W(T))- 2 + 2max|A i pao+)| < 

i 

Kin(W(T))- 2 + 4,X^{W{T))- l ^/bK a , b K D , K s{I). 

On the right-hand side only K a ^ and \ m m{W(T))~ 2 depend 
on the temperature, both are decreasing in T and go to zero as 
T goes to infinity, i.e., 

gives an implicit equation for the temperature T c above which 
there is no bi-partite entanglement left in the system. 

VII. NEAREST-NEIGHBOR INTERACTIONS 

In this section we consider nearest-neighbor interactions 
and periodic boundary conditions in D spatial dimensions, 
i.e., block-circulant matrices V. For n — > oo this is a special 
case of the matrices considered in Section[V]as boundary con- 
ditions become irrelevant in this limit, i.e., the upper bound 
coincides with the one derived for finite-ranged interactions. 
For a tighter bound see Appendix 151 We will now make use 
of the circulant structure of V to show that for matrices of this 
kind it is possible to also derive a lower bound on the entropy 
of entanglement that it is proportional to the surface area of 
/. We write M = circ(M) for the circulant matrix M whose 
first block column is specified by the tupel of matrices M. We 
can then recursively define V = Vn via 

Vs+i = circ(Vs,-cl n s,O n s,. . . , n s , -cl n «), 

where Vi = circ(l, — c, 0, . . . , 0, — c). We choose an energy 
scale in which Vi.i = 1 and we demand c > l43l . This 
circulant structure leads to the following properties of V: the 
eigenvalues of V are given by 

D 

Xi(V) = 1 - 2c cos (2-Kis/n) , i e C, 

6=1 

in particular the maximum eigenvalue is given by X mllx (V) < 
1 + 2cD, where equality holds for n even. The minimum 
eigenvalue reads A m i n (V) = 1 — 2cD, i.e., positivity of V de- 
mands c < 1/2D. Note that the assumption that c < 1/2D 



independent of / and O is essential for the argument. If we 
allow for an n-dependence of c as it arises, e.g., in the field- 
limit wherec = l/(l/n 2 + 2_D), then in one spatial dimension 
one will encounter an area law up to a logarithmically diver- 
gent correction. This behavior - resembling the behavior of 
critical spin chains and quadratic fermionic models - will be 
studied in more detail elsewhere. 

Furthermore the circulant nature of V yields the following 
explicit expressions for the entries of V ±x l 2 



V 



±1/2 



^-rrifcfi— j) J n 



(Xk(V)) 



±1/2 



(14) 



This reasoning leads to the following properties for i ^ j that 
are crucial for the present proof (see AppendixlBl. 



[V 



-l/2l 



> 



2d(i,j) 



d(i,j) 



(2d(i,j)-l)Uti\ i s-h\ 



Substituting these results into the general lower bound dlOt 
and keeping only terms with d(i, j) = 1 immediately yields 



E 



where 



log[V^]c 2 



l + 2c£> 
1 - 2cD 



(15) 



in this case. This generalizes the result of Ref. 0] to regions 
/ of arbitrary shape. We expect that these lower bounds can 
be generalized to other interactions V and numerical results 
suggest that these bounds hold quite generally. 



VIII. SQUARED INTERACTIONS 

A simple special case is related to a certain kind of inter- 
action: this the one where the interaction matrix V can be 
written as 

V = M 2 

with a real symmetric matrix M corresponding to a finite- 
range interaction. Then, the covariance matrix associated with 
the interior is nothing but 



7o 



M~ 



© M\j. 



(16) 



In this case, the symplectic spectrum can be determined in a 
fairly straightforward manner. We mention this case also as 
it appears to be an appropriate toy model as a starting point 
for studies aiming at assessing the symplectic spectrum of the 
reduction itself and therefore the spectrum of the reduced den- 
sity matrix l26ll . 
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We are looking for the spectrum of the v(I) x v(I) matrix 
Aj as in Eq. © which now takes the form (recall that we label 
the entries of A] by vectors i,j £ /) 



'1,3 



R 



h3 ■ 



k£0 



It is now the central observation that the number of rows of 
the matrix R that are nonzero - and therefore also the number 
of eigenvalues that are nonzero - is proportional to the surface 
area of /, as M is a banded matrix. That is, 

Mij = 

for d(i,j) > k/2, meaning that Ri j = if d(k,j) > k/2 
for all k £ O. As the eigenvalues of Aj are bounded through 
the pinching inequality by 1 and A max (W / )/A m i n (M / ), an up- 
per bound linear in the surface area of the interior follows im- 
mediately from the fact that the number of eigenvalues enter- 
ing the sum in Eq. is proportional to s(I). This argument 
demonstrates that in this simple case, one immediately arrives 
at bounds that are linear in the number of contact points. 

The task of finding the eigenvalues of Aj explicitly is now 
reduced to finding the eigenvalues of the sparse matrix R, 



Hi(Ai) = ^l-Xi(R). 

This case is particularly transparent in the one-dimensional 
case, D = 1, and for M being a circulant matrix with first 
row (1, — c, 0, . . . , — c) for < c < 1/2. The potential ma- 
trix V corresponds then to nearest-neighbor interactions, to- 
gether with next-to-nearest-neighbor interactions. In this one- 
dimensional setting, we set / = [1, . . . , to]. The matrix R 
then takes the simple form (cf. Fig.|4} 




FIG. 4: Matrix entries of R for the case D = 1, C = [1, 100], 
/ = [1, 15], V = M 2 , where M is a nearest-neighbor circulant 
matrix. Note that for this particular form of V the entries Rij are 
exactly zero for j ^ 1, m. This is in contrast to the general form of 
R depicted in Fig. [2] The inset shows the entries Tl^i as bars and 
log(i?i,i) in blue (cf. Eq. {T7J). 



From the circulant structure of M we have 

I JL^ e 2mk{i-j)/n 

n 1 - 
fe=l 



[m- 1 ; 



1 - 2ccos(27r£:/n) 
1 x--^ cos(27rfc(i — j)/n) 

k=l 



1 x - CC 

n^r-2ccos(27rfc/n)' 



k—m-\-l 



-c ([M-\ m+1 6 jim + [M-^.n^i) 



To find its eigenvalues, we calculate det(i? — Al), which is 
straightforward for a matrix of this form. 



det(i?- Al) = -(-A) 



m-2 



dct 



Rll — A Rlr 



R 



71 1 1 



A 



— (— A) m 2 (Ri m R m i — (Rn — X)(R mm — A)), 

reflecting the fact that m — 2 eigenvalues are zero, i.e, the 
number of nonzero eigenvalues is s(J) = s([l, . . . , to]) = 2. 
We find for the non-vanishing eigenvalues 



2A± — Rn + Rmm 

± ( (Rll — Rmm) 2 + 4i?i m i? m i 



1/2 



Symmetry of M 1 yields 

A± = Ru ± Rim = 



([M-\ n ± [M~\ 



i.e., for large n we arrive at 



[M~ 1 ]i n 



1 



1 



2cVl - 4c 2 2c' 

1 ( 1- \/l - 4c 2 ^ 



2c 



(17) 



Note that these expressions are asymptotically independent of 
to as shown in the inset of Fig.|4] In this limit we finally arrive 
at 



A 



1 



1 



2 2V1 - 4c 2 



(18) 



This expression specifies the symplectic spectrum of the re- 
duction in a closed form. 



IX. DIVERGING CORRELATION LENGTH AND 
AREA-LAW OF THE ENTROPY 

In this section we will analytically demonstrate that there 
exist Hamiltonians for which the ground state two-point cor- 
relation functions diverge whereas the geometric entropy is 
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still bounded by an area-law. In spin systems with a long- 
range Ising interaction such a behavior has been observed in 
the one-dimensional case and sketched for higher dimensions 
I23I1 . Here, we present a class of examples of this type valid 
in arbitrary dimension D, for which one can prove the validity 
of an upper bound linear in the boundary area. Moreover, this 
set of examples is not restricted to cubic regions. The interac- 
tion is here a suitable harmonic long-range interaction. This 
analysis shows that a divergent two-point correlation function 
alone is no criterion for a saturating block entropy in the one- 
dimensional case, and for an area-dependence in higher di- 
mensions. 

Consider the matrix M with entries 



Mi 



1 



d(i,j)<*' 



for some a > 0. Now set V = A/~ 2 . This choice implies that 
the correlation function Gi j = [V -1 / 2 ]^ decays only alge- 
braically. We will show that despite this fact one still has an 
upper bound linear in the boundary area of / for appropriate 
values of a. 

Firstly, we have to make sure that the maximum eigenvalue 
of V can be bounded from above independent of n. From 
Gershgorin's theorem (see, e.g., Ref. 13910 we know that for 
every eigenvalue X(M) there exists an i such that 



\X(M)-Mi,i\ <J2 M iJ> 



jec 



smaller values of a, such a behavior can be expected. Steps 
towards tightening these bounds seem particularly feasible in 
case of cubic regions, where we conjecture that for a > D 
we arrive at an area-dependence. This analysis shows that for 
long-range interactions, an area-law in the degree of entangle- 
ment can be concomitant with divergent two-point correlation 
functions. 



X. AN AREA LAW FOR CLASSICAL CORRELATIONS 

In previous sections we have considered the entanglement 
between some region and the rest of a lattice of interacting 
quantum harmonic oscillators with Hamiltonians that are at 
most quadratic in position and momentum. We have demon- 
strated that both for the ground state and the thermal state of 
the entire lattice the quantum correlations, i.e., the amount of 
entanglement, between the region and the remainder of the 
lattice is bounded by quantities proportional to their boundary 
surface area. 

This suggests similar questions concerning the classical 
correlations between a region and the rest of the lattice in the 
corresponding classical systems when the lattice system as a 
whole is prepared in a thermal state. In this section we will 
demonstrate that indeed analogous area laws hold. We will 
note furthermore that there is a quite striking intimate rela- 
tion between the classical system with potential matrix V c and 
the quantum system with a potential matrix V q = V 2 in this 
context. 



A. Hamiltonian, entropy, and mutual information 



i.e.. 



Amin(M) > min 



Mi 



5>, 

dec 



and therefore 

b = A max (y) = A max (A/- 2 ) 



< 1. 



Substituting this and the specific form of M into the general 
expression for the entropy, we obtain 



ln(2) 



E 



1 00 1 

'? - E uu AT 5 = Et^ 



1=1 



< 2 2D s(I)J2l 2D ~ a ~ 1 , 



1=1 



which converges for a > 2D to 

o2D+2 



Ef < 



ln(2) 



■C(a-2£> +!>(/), 



(19) 



where ((x) is the Riemann zeta function. Note that these 
bounds are not necessarily tight in the sense that even for 



For the following considerations we use the classical equiv- 
alent of the quantum mechanical Hamiltonian Q, namely 
again 



3 3 



1-3 



(20) 



where now x = [x\, . . . , x n n) and p — (pi, . . . ,p n r>) are 
vectors of classical position and momentum variables, respec- 
tively, and V denotes again the potential matrix. The state of 
the classical system is characterized by a phase space den- 
sity g = £>(£), a classical probability distribution, where 
£ = (xi, . . . , x n D , pi , . . . , p n D ) denotes all the canonical co- 
ordinates in phase space. For nonzero temperatures T > 
this phase space distribution is given by the Boltzmann distri- 
bution 



9(0 



1 



-mm 



(21) 



where (3 — 1/T (as before, k = 1), and 



Z := / d£e 



-em) 



Given this density in phase space we will encounter a famil- 
iar ambiguity when defining the entropy of the system. Using 
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the discrete classical entropy 14311 - we split the phase space 
into cubic cells each with a volume h 2N , where N = k D and 
h > being an arbitrary constant. From the phase space den- 
sity we obtain the probability associated with each of these 
cells which in turn can be used to determine the entropy func- 
tion of this probability distribution. We will now make use of 
the multiple indices % = (ii, . . . , Izv) and j = (ji, . . . , jjy), 
assuming that the cell corresponding to is centered 

around x = (hi%, . . . , Uin) and p = (hjx, ■ ■ ■ , hjN). That 
is, for each degree of the n D degrees of freedom, the phase 
space is discretized. The contribution of each cell to the dis- 
cretized probability density is then given by 



P(i,j)= I dZQ{t-) 

J Cell 



As usual, the discrete classical entropy is then defined as the 
corresponding Shannon entropy 



(22) 



We will denote the discrete classical entropy with respect to 
the degrees of freedom of the interior with Si(h), the entropy 
of the exterior with So (h). 

The value for the entropy will depend on the choice of h 
and in the limit h—*0 this entropy definition will diverge due 
to a term proportional to — \og(h 2N ). In classical statistical 
mechanics this problem is avoided with the help of the third 
law of thermodynamics. The entropy itself is however not 
the quantity that we wish to compute but rather the mutual 
information Ij between the interior I and the rest of a lattice, 
denoted as before by O. This classical mutual information 
meaningfully quantifies the classical correlations between the 
inner and the outer. In that case we find that the limit h — ► 
exists and that the mutual information I can be defined as 

h := lim (Si{h) + S Q (h) - S c (h)) . 

Following these preparations we are now in a position to de- 
termine the mutual information between a region and the rest 
of the lattice explicitly when the lattice as a whole is in a ther- 
mal state. 



B. Evaluation of the mutual information 

For the evaluation of the mutual information we need to 
determine the entropy of the total lattice C, as well as the 
entropy determined by the reduced densities describing the 
two regions I and O. To this end we carry out the partial 
summation over all degrees of freedom of region O in order 
to find the reduced phase space density qj describing region 
/ only. Employing the Schur complement we find that the 
reduced density gj is described by the Boltzmann distribution 
corresponding to the same temperature and the Hamiltonian 



H, 



iei 



(23) 



An analogous result holds for the reduced phase space density 
of region O. 

For a thermal phase space distribution Eq. d2 1 1 correspond- 
ing to a classical Hamiltonian function of the form Eq. (I20> 
we can compute the entropy straightforwardly, to find 



S A = -ilogdet(V| A ) 1 +v(A)\og 



2tt 
J 



v(A) 



A G {/, O, C}, which increases with temperature as ex- 
pected. For the mutual information we find 



1 det V Ldet V L 
- log - — . 

2 B det V 



which is, perhaps surprisingly, independent of temperature. 
Using Jacobi's determinant identity 



det V\ detV~ 1 = det V^ 1 



this expression can be rewritten as 



I i = \ lo S det 



v\ r v- 



= -log det [l-.fi] 



It is now advantageous to notice the close connection of this 
expression, in particular of the matrix R 



R 



hi 



keo 



with those that arise in the quantum mechanical problem that 
we have treated previously in Section IVIIII Indeed, the clas- 
sical problem for a system with potential matrix V c is related 
to the quantum mechanical system with the squared poten- 
tial matrix Vq = V 2 . This formal similarity arises because in 
Section rvilll we have shown that for a lattice of quantum har- 
monic oscillators with potential matrix V q = V 2 in its ground 
state the symplectic eigenvalues of the covariance matrix de- 
scribing region I alone are exactly the standard eigenvalues of 
the matrix (1 — i?) 1 / 2 . The properties of these eigenvalues 
have already been discussed in detail in Section IVIIII This 
allows us now straightforwardly to establish the area theorem 
for the mutual information in a classical system employing the 
result for the corresponding quantum system. 

This establishes in particular that the classical correlations 
as measured by the mutual information Ic,i between the dis- 
tinguished region / and the rest of the lattice O = C \ I satisfy 

<k ■ s(I) < Ii(T) < 4 • s(I) 

for large v(I) and appropriate constants c[ 2 > independent 
of O and /. (Compare also the assessment of the thermo- 
dynamical entropy of parts of classical fluids in Ref. 14411 .') 
In summary we have seen that the area-dependence of corre- 
lations is not restricted to quantum systems, as long as one 
replaces the notion of entanglement - representing quantum 
correlations - by the notion of classical correlations in a clas- 
sical system. 
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XI. SUMMARY AND OUTLOOK 

In this paper, we have considered the question of the area- 
dependence of the geometric entropy and the distillable entan- 
glement in general bosonic harmonic lattice systems of arbi- 
trary dimension. The question was the general scaling behav- 
ior of these measures of entanglement with the size of a dis- 
tinguished region of a lattice. Such an analysis generalizes as- 
sessments of block entropies in the one-dimensional case. Us- 
ing methods from entanglement theory, we established bounds 
that allow for a conclusion that may be expressed in a nutshell 
as: in surprising generality, we find that the degree of entan- 
glement scales at most linearly in the boundary area of the 
distinguished region. This analysis shows that the intuition 
that both the interior and the exterior can be approximately 
disentangled up to a layer of the thickness of the two-point 
correlation length by appropriate local unitaries carries quite 
far indeed. 

For cubic regions / = [1, ...,m} xD the area law can be 
formulated as 



E? 



(24) 



where O is the Landau theta. 

Such area-laws are expected to have an immediate implica- 
tion on the accuracy to which ground states can be approxi- 
mated with matrix-product states and higher dimensional ana- 
logues in classical simulations of the ground states of quantum 
many -body systems 11811 . After all, the failure of DMRG algo- 
rithms close to critical points can be related to the logarithmic 
divergence of the block entropy in the one-dimensional case. 

The findings of the present paper raise a number of interest- 
ing questions. Notably, in general quantum many-body sys- 
tems on a lattice (fermionic or bosonic), what are necessary 
and sufficient conditions for an area-law in the above sense 
to hold? Clearly, as we have seen above, the divergence of 
two-point correlation functions alone is not in one-to-one cor- 
respondence with an area law. It would be interesting to con- 
sider and possibly decide the conjecture that a one-to-one re- 
lationship between a system being critical and not satisfying 
a law of the form as in Eq. i24i holds if one (i) restricts atten- 
tion to systems in arbitrary dimension with nearest-neighbor 
interactions, and (ii) grasps criticality in terms of two-point 
correlation functions with algebraic decay, concomitant with 
a vanishing energy gap. Note that the latter two criteria of crit- 
icality have not been rigorously related to each other yet and 
may indeed not be simultaneously satisfied in lattice systems. 
The general relationship is still awaiting rigorous clarification. 

As steps towards such an understanding of a relationship 
between criticality and properties of ground state entangle- 
ment in more than one-dimensional fermionic and bosonic 
systems, it seems very interesting to study models different 
from the ones considered in this paper. For example, comple- 
menting our bosonic analysis, area-laws in fermionic critical 
systems have been addressed l24H25ll . where logarithmic cor- 
rections have been found also in higher-dimensional settings. 
Other settings in the bosonic case, corresponding to field the- 
ories beyond this quasi-free setting, are also still not clarified. 
In particular, the scaling behavior of the geometric entropy 



in general bosonic theories in higher dimensions is far from 
clear. Then, the case of finite-size effects in harmonic lattice 
systems where the correlation length is larger than the full sys- 
tem, resembling the critical case, will be presented elsewhere. 
In less generality, it seems also feasible to identify the prefac- 
tors of the leading and next-to-leading order terms in an area- 
law of the geometric entropy. It is the hope that the present 
work can contribute to an understanding of genuine quantum 
correlations in quantum many-body systems and inspire such 
further considerations. 



XII. ACKNOWLEDGMENTS 

We would like to thank C. H. Bennett, H. J. Briegel, 
J. I. Cirac, W. Diir, M. Fleischhauer, B. Reznik, T. Rudolph, 
N. Schuch, R. G. Unanyan, R. F. Werner, and M. M. Wolf 
for discussions. This work has been supported by the EP- 
SRC QIP-IRC (GR/S82176/0), the EU Thematic Network 
QUPRODIS (IST-2002-38877), the DFG (SPP 1078 and SPP 
1 1 16), The Leverhulme Trust (F/07 058/U) and the European 
Research Councils (EURYI). 



APPENDIX A: EXPONENTIAL DECAY OF ENTRIES OF 
MATRIX FUNCTIONS 



The result concerning the exponential decay of entries of 
matrix functions of Ref. Bill relies on the fact that the p-th 
power A p of a fc-banded matrix A = (Aij) is pfc-banded, i.e.. 



for \i — j\ > pk/2 for a matrix A with Aij = 



[AP] 

for \i — j\ > fc/2. For the purposes of the pres ent p aper, we 
will need a generalization of the result of Ref. 14111 to block 
banded matrices. We refer to V — (Vi.j) as being fc-banded, 
if Vi.j = for d(i,j) > fc/2. It can be proven by induction 
over p that the p-th power of V is pk banded in this sense. 
This enables us to formulate the general form of Ref. EUl as 
follows. 



Let V 



(Vi, 



be a fc-banded symmetric matrix, i.e., 
Vij = for J2s N<5 ~ Js\ > Define a = X min (V), 

b = A raax (y), k = b/a, 



(b — a)z + a + b 



and e x as an ellipse with foci in —1 and 1 and half axes a, (3, 
X = a + (3. 

Now let / : C — ► © be such that fotp is analytic in the interior 
of the ellipse e x , x > 1- an d continuous on e x . Furthermore 
suppose (foip)(z) € IRforz £ R. Then there exist constants 
K and q, < K, < q < 1 such that 



[f(V)]ij <Kq^^-^K 
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where 



K = max < max Aj(/(V)) 



2/fc 



2xM(x) | 
X-l J 



3= - , M( X )=max|(/o^)(z)|. 

To bound the entries of u>~, we apply the above theorem to 
the function 



u-{z) 



3 v^/T _ 1 



For 1 < x < (y 7 ^ + 1) 2 /( K ^ 1), we have that o ?/> 
is analytic in the interior of the ellipse e x and continuous on 
e x , i.e., satisfy the assumptions of the above theorem. To 
make K and q specific, we choose 



(^+1) 2 K+l 1 



K — 1 



K- 1 



which yields 



where 



k — 1 

K+l 



2/ ft 



e^ T - 1 k + 1 
e^/ 7 + 1 i] 



V 



K 



a- 



1/2 



and for zero temperature K a _b reduces to 

K+l 



K a ,, 



V 



These findings explicitly relate the spectral properties of the 
Hamiltonian to the two-point correlation functions. 



APPENDIX B: ENTRIES OF THE CORRELATION MATRIX 
FOR THE NEAREST-NEIGHBOR CASE 

In this appendix, we make the evaluation of the entries of 
V 1 ! 2 and I/ -1 / 2 of the important case of nearest-neighbor in- 
teractions specific. The power series expansion of the square 
root is given by 



k=l 



a k =11-21-' 

i=i 



H ~ 2k -1 ~ a « ~ 2*(2fc-l)' 

which is valid for |x| < 1, i.e., the positivity constraint < 
c < 1 /2D allows us to write 



(A fe (y)) ± * = l + ^a ; ± (2c^cos( 



1=1 x s=i 

D 



1=1 M=l 



Using the multinomial theorem, we have 



5=1 



2^ _2 ll k jL \ns 

D I e " + e 



5=1 



ngl 



where the sum runs over all rig with J2s ns = Now, apply- 
ing the binomial theorem 



Tig \ 2xifca(n a -2r) 

e 

r 



( e ^+ e "^J =E 

r=o 

Substituting all the above into Eq. (1 1 41 . we find 



1=1 



6=1 



[V ±1/2 ]ij = 8ij + £ 4c l E 11 IT f(ds,m), 
where dg = ig — js and 



f(ds,n s ) = ^2 



E 



- 1 27ri/c(rf ( 5+n ( 5 — 2r)/n 



— ' (nx — r)!r! 

r=0 v ; fc=l 



To the sum over r only terms with ds + ns —2r = zn for some 
z E % contribute. We thus arrive at the following expression 
for the entries of V ±x l 2 

oo D n§ 1 

1 r=0 V d ' 



1=1 Es n «= ! 



i.e., for nearest-neighbor interactions the product 
— \V~ 1 / 2 ]ij\V +1 / 2 ]i.j is always positive for i =/= j. 
This does not hold in general and makes it difficult to obtain 
explicit bounds on the entries of V ±x / 2 for more general 
interactions V. To obtain a lower bound we keep only the 
term I = d(i,j) and rig = \dg\. The restriction on r is then 
fulfilled for r = (r = \d s |) if d 5 < (d s > 0), yielding 



\[v ±1/2 \i,\ > *%jf d{id) 



d(i,j) 



Us=i \n-js\ 



d(i,j) 



(2d(i,j)-l)Uti\^-3s\ 



It is also possible to obtain an upper bound that is tighter than 
the one derived for finite-ranged interactions: The elements of 
y± x / 2 are symmetric under i — j — > ne$ — (i — j), where 
eg is a unit vector along dimension 5. Thus, we can demand 
—n/2 < is — js < n/2. Then we find that ns has to be larger 
or equal to \ig — jg\ otherwise the restriction on r can not be 
fulfilled. This in turn means that I has to be larger or equal to 
d(i, j). We then obtain an upper bound by summing all terms 
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FIG. 5: Visualization of the enumeration of Ni as in Appendixlcl As 
before, oscillators belonging to the distinguished (shaded) region I 
are marked •, the outside ones are shown as O. dO is shown as the 
orange shaded area. Ma(%) = {j £ C\d(i,j) < 4} for a certain 
oscillator • is shaded green, its surface oscillators 7714(2) = {j £ 
C\d(i,j) = 4} are depicted by •. 



in the sum over r regardless of the given restriction, yielding 

D 



i[^ ±1/2 k,-i < E 



l>d(i,j) £ 5 «<5=i S=l 

e a ^ 2c °y 
i>d(ij) 



= E«t, (i ,i)( 2 ^) i+dW) 

1=0 

00 

< 4 (i>j) (2cD) d ^J2( 2cD ) 1 



1=0 



a d(i,3) 

1 - 2cD 



(2cD) 



d(i,i) 



APPENDIX C: ENUMERATING THE RELEVANT TERMS 
IN THE AREA LAW 



We start by identifying the set of oscillators that can con- 
tribute to Nu I > 1, 



w = E E L 



j'eo iei 
d(i,j)=i 



Oscillators j e O can only contribute if their distance to the 
boundary dO is not larger than I — 1, 



dO = {jeO\3ieI:d(i,j) = l}. 
Thus, we can restrict the sum over O to the set A\ 



A l = U {3 £0\d(i,j)<l-l}, 

iedO 

i.e., we can write 



*'=E E ^E E i<i^k, 

oeAi iei oeAi iec 

d(i,o)=l d{i,o)=l 



where m; is the number of surface oscillators of a ball with 
radius / within the metric d, i.e., mj < 2(21 + l) ^ 1 for 
I > 1. Using the fact that \dO\ < Ni — s(I), \A t \ can now 
be bounded from above in the following way 

\A t \ < \dO\Mt_x < s(J)M,_i, 

where M; is the volume of a ball with radius I within the met- 
ric d, i.e., Mi < (21 + 1) D . To summarize, we have 



Ni < 2(21 - l) u (2l + l)^ 1 ^/) < 2(2l)* 1J - L s(I) 
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